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Abstract. We present the next-to-leading strong interaction corrections to the AF = 2 Hamilto- 
nian in the MSSM with exact diagonalization of the squark mass matrices. These results allow 
phenomenological studies of neutral meson mixing in scenarios with non-degenerate squarks, with 
control over the renormalization scale and scheme dependence. 
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Flavor changing low-energy processes such as weak meson decays or the mixing of 
neutral mesons, are unique for indirect new physics searches because they are very sen- 
sitive to heavy degrees of freedom in loop amplitudes. This is due to the suppression of 
flavor-changing neutral currents in the SM, in contrast to the enormous flavor violation 
typical of new physics scenarios with generic flavor structure. Precise measurements - 
specially those related to K and B mesons- have put stringent constraints on either the 
new physics scale or the flavor non-universality of the new couplings [1]. This is explic- 
itly clear in SUSY models, and in particular in the MSSM, which is the paradigmatic 
extension of the SM. 

In the MSSM, the leading source of flavor violation arises from a misalignment be- 
tween quark and squark mass eigenbases induced by soft SUSY-breaking terms: in the 
basis for the superfields in which the quark mass matrices are diagonal (the so called 
super-CKM basis), the squark mass matrices contain non-zero non-diagonal elements 
(mass insertions), which are responsible for the flavor transition. Extensive phenomeno- 
logical studies have been performed which constrain strongly the size of these mass 
insertions, based on FCNC processes 10, SB SB], vacuum stability requirements iSB] 
and charged-current processes [9]. 

However, the increasing precision with which these quantities must be determined 
requires theoretical calculations of flavor violating processes at increasing order in per- 
turbation theory, as well as more precise determinations of non-perturbative parameters. 
Here we focus on the calculation of matching conditions (Wilson coefficients at the 
SUSY scale) for the AF = 2 Hamiltonian, which are required for the theoretical deter- 
mination of K — K, D — D and ^ ^ mixing amplitudes. The case of Bs—Bs mixing 
has recently caused some excitement in relation with a deviation of the measured mixing 
phase relative to the SM expectation [n^. [lll [l2]. and if such deviation is confirmed, a 
precise calculation of that phase in different models will be very convenient. 

Leading order (LO) strong interaction matching conditions in the MSSM have been 
computed in Refs. lfl3l [3, [isl], and arise from squark-gluino box diagrams (Fig. [IJa). 
The corresponding next-to-leading order (NLO) corrections arise from two loop dia- 
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FIGURE 1. A sample of the diagrams contributing to the MSSM matching conditions for — 
mixing at (a) leading order, and (b) next- to-leading order in 



grams such as those shown in Fig. [TJb, and have been computed in Ref. [16] within the 
degenerate Mass Insertion Approximation (that is, at the leading order in an expansion 
in mass insertions, and assuming that the diagonal elements of the squark mass matrix 
are equal). Here we report on the calculation of the NLO strong interaction corrections 
for non-degenerate squark masses, that is, with exact diagonalization of the squark mass 



matrices. This calculation was performed in Ref. [ITQ. extending the known results be- 
yond the Mass Insertion Approximation (MIA), and the full expressions for the Wilson 
coefficients can be found in that paper. 

These NLO corrections valid for non-degenerate squark masses are important for two 
reasons: 

1. At LO, the renormalization scale and scheme cannot be specified for the strong 
coupling, leading to a large uncertainty related to scale and scheme ambiguities. At 
NLO, however, scheme independent results and NLO scale invariance are accom- 
plished. The reduction of this uncertainty is numerically important due to the large 
anomalous dimensions of the AF = 2 operators involved. The analysis of Ref. [,161 
shows that this uncertainty is reduced from 10-15% down to a few percent. 

2. Scenarios with large squark mass splittings can give rise to a quite different pat- 
tern of correlations between AF = 1 and AF = 2 observables, such that AF = 2 
processes can be affected very differently by the experimental bounds on AF = 1 
processes. An example of how this comes about in a hierarchical scenario is given 
in Ref. [fisll , where it is shown that the bounds from 5 — > Z^y can be partially evaded 
allowing a large Bs mixing phase. 

From the results with exact diagonalization of the squark mass matrices presented 
here and in Ref. it is possible to recover the Wilson coefficients in the MIA. This is 
done by relating the rotation matrices to the mass insertions, expanding to lowest order in 
the mass insertion expansion and taking the limit of equal diagonal elements of the mass 
matrices. The corresponding expressions agree with the results of Ref. il<^ . verifying 
the correctness of their calculation. Moreover, in order to provide a clear comparison 
between the degenerate and non-degenerate scenarios, it is convenient to expand the 
exact results to lowest order in the mass insertion expansion but keeping the squark 
masses non-degenerate. This scheme might be called the non-degenerate mass insertion 
approximation (NDMIA). 

In order to get an idea of the size of the new contributions, we consider a particular 





FIGURE 2. Relative importance of the full NLO result with respect to the LO, for the Wilson coef- 
ficients C2 and C3, as a function of = rh\2/mg, where mi2 is a common mass for first and second 
generation squarks. The dashed lines correspond to the degenerate (MIA) scenario. We have chosen 
400 GeV for third generation squark masses and nig = fi — 350 GeV, where /i is the matching scale. 



scenario in the NDMIA for the case of — mixing. This scenario is based on the 
"hierarchical" setup of Ref. lIlSll . in which the first two generation squarks are given a 
common mass mi2, different from a common mass for the third generation squarks. 
Moreover, is assumed to be near the electroweak scale in order to ensure naturalness, 
while mi2 could be well above (thus reducing the conflict with kaon physics data, for 
example). 

The plots in Fig. |2] illustrate the relative importance of the NLO corrections, as a 
function of the mass splitting between the light and heavy squarks. The NLO correction 
is typically a ~ 10% effect in the degenerate case, but its importance increases with the 
mass splitting. For heavy squarks of about a TeV, the NLO contribution to C3 can be up 
to a ~ 25% correction. 

In order to analyze more closely the role of the mass splittings in the NLO corrections, 
we consider C\,C^ and C5 as a function xi = m^/m^, for different splittings between xi 

and Xh = mlj/tn^- This is shown in Fig.[3l where the dashed lines correspond to xi = Xh 
(that is, the degenerate case), and -departing smoothly from that limit- the solid lines 
show increasing values of x^/xi. In these plots we take the dimensionless mass insertions 
Sl,r = Srx = and nig = 11 = 350 GeV. We see that increasing the heavy scale tends to 
reduce systematically the size of the NLO contribution, being largest in the degenerate 
case. 

However, caution must be taken when interpreting these results, since the NLO Wil- 
son coefficients are renormalization scale and scheme dependent. These results corre- 
spond to the renormalization scheme adopted in Ref. [17] (NDR with some subtleties), 
and will vary in other schemes. Of course, the scheme dependence cancels in physical 
amplitudes (against the scheme dependence of the matrix elements of the operators), 
so the true impact of these NLO corrections can only be established by analyzing their 
effect on observables. A full phenomenological analysis of these corrections and their 
impact on the bounds on the mass insertions (beyond the mass insertion approximation) 
will be presented in the future. 
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FIGURE 3. NLO Wilson coefficients Cf^^, C^^^ and C^^^ as a function of x/ = m^m?, where is a 
common mass for third generation squarks. The plots are in units of {a^ / T^)8l^ and {oc^ / 7:)5ll5i{i{ for Ci 
and C4,5 respectively. The different lines correspond to different values of = mj2/'"|' ™d range from 
the MIA case, X/, =x/ (dashed) tox/, = 1.5x/,2x/,3x;,4x/,5x/,6x/. 
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